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^ \ Abstract 

^vq ' We investigate tlie observable consequences of creating cosmic texture by 

f^ \ breaking a global SU(3) symmetry, rather than the SU(2) case which is generally 

t^^ ' studied. To this end, we study the nonlinear sigma model for a totally broken 

^1* \ SU(3) symmetry, and develop a technique for numerically solving the classical 

r^ ' field equations. This technique is applied in a cosmological context: the energy 

Qh! of the collapsing SU(3) texture field is used as a gravitational source for the 

production of perturbations in the primordial fluids of the early universe. From 
q_) \ these calculations, we make predictions about the appearance of the anisotropics 

i-S^ ' in the cosmic microwave background radiation (CMBR) which would be present 

if the large scale structure of the universe was gravitationally seeded by the 
collapse of SU(3) textures. 



1 Introduction 

Over the last few years, there has been a great deal of research on the effects 
of topological defects in the the early universe. In particular, one class of topo- 
logical defect, textures, has emerged as an interesting canidate for creating the 
seeds of large scale structure in the early universe Q . In three spatial dimen- 
sions, textures can be created in any symmetry breaking scheme where the 
vacuum manifold has a nontrivial third homotopy group, -k^. This situation is 
quite generic - for example if any nonabelian Lie group is completely broken, 
773 is nontrivial. There are a multitude of symmetry breaking schemes which 
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could produce textures (a general classification has been given in |^). But so 
far the only case which has been studied in any detail has been the simplest 
case, where the vacuum manifold is SU(2) ~ S^. An accurate nonlinear sigma 
model code has been developed for SU(2) textures al.^, and it has been pos- 
sible to extract detailed predictions of the power spectrum of inhomogeneities 
and nongaussianity of the CMBR anisotropics , [ i| . 

In the light of such investigations, it is natural to ask whether the primary 
results from these simulations are particular to SU(2), or whether they are more 
general, pertaining to many or all field theories which would create texture. To 
answer this question, one must study the effects of textures that come from at 
least one other group. All simple Lie groups have a nontrivial tts so in principle 
one might choose any one of them. However, simplicity dictates SU(3) as the 
obvious choice. In this paper we shall study the effects of textures coming from 
a broken global SU(3) symmetry upon the CMBR anisotropics and density fluc- 
tuations of the early universe. A complementary investigation of more general 
texture theories has been recently conducted by Sornberger et. al. [||. 

There is also some motivation from particle physics to study the effects of 
a totally broken global SU(3) symmetry. Joyce and Turok, |^, have pointed 
out that if the family symmetry among the quarks and leptons is broken by the 
Higgs mechanism, then we must have Higgs fields which totally break some three 
dimensional representation of a simple Lie Group. The only simple Lie groups 
which have three dimensional representations are SU(2) and SU(3). Joyce and 
Turok argue that SU(2) does not work well, as it does not easily break to make 
one quark family much more massive than the others, whereas in SU(3) this 
occurs naturally. A gauged SU(3) family symmetry would be anomolous, so such 
an SU(3) family symmetry must be global. When this symmetry is broken, the 
Goldstone modes of the Higgs fields would produce cosmic texture. In particle 
physics terms, the point of the present paper is to explore the cosmological 
consequences of a global SU(3) family symmetry broken at the GUT scale. 

In this paper, a technique for studying textures from a broken global SU(3) 
symmetry is developed and applied, with the goal of making predictions about 
fluctuations in the cosmic microwave background. In section y, we state the 
problem and derive the relevant nonlinear sigma model equations y]. In sec- 
tion y, we develop a finite difference algorithm based upon these equations, and 
in section we describe tests of the stability and accuracy of this algorithm. 
In section IS we present some results about the scaling behavior of this algo- 
rithm, and the appearance of the CMBR fluctuations if the primordial density 
fluctuations of the universe come from this type of texture. 

We find that the cosmological consequences of textures coming from a totally 
broken SU(3) symmetry are in most respects very similar to those from SU(2) 
textures. We present some preliminary evidence of differences in the shape of 
the CMBR anisotropy power spectrum, but more accurate techniques will be 
needed to fully resolve these. 



2 Formulation of the problem 

We wish to totally break the fundamental representation of SU(3) via the Higgs 
mechanism. In M, Joyce and Turok describe this process in detail. In order 
to do this, our Higgs field must consist of two complex triplets, S^,S^ £ C^, 
which have a "Mexican hat" potential whose minimum is topologically equiv- 
alent to SU(3). Without loss of generality, one can choose these fields so that 
the minimum of the potential is described by: 

i:^^-E^ = vl ^^^i:^^vl si^s^^o. (1) 

Here vi , f 2 are the vacuum expectation values of S , S , and are not necessarily 
equal to one another. (See equation (62) of |^.) Any particular set of S , S 
satisfying (|l|) clearly define a unique member of SU(3), for if we define 

^ = — , X-—, (2) 

then we can immediately construct an SU(3) group element, in the fundamental 
representation, 

U^{ iP^xx^ tP X \ ■ (3) 



Note that the vacuum manifold is only identical to SU(3) geometrically if 
Vi = V2- Otherwise it represents a 'stretched' version of SU(3). 

We wish to construct a nonlinear sigma model which will describe the classi- 
cal dynamics of these fields, assuming that they are restricted to this minimum. 

The first thing we need are the classical field equations. The Lagrangian 
density for this system is 

C = vlWf,xp^W''ip + vlWf,x^W''x + 

A^(V' V - 1) + ^xix^X - 1) + Av^x^^X + A;^xV , (4) 

where the As are all Lagrange multipliers, with A^,Aj^ real, and A^^ complex 
valued. Note that the theory has one free parameter, the ratio of VEVs of 
the two Higgs fields, ui/w2- Using this Lagrangian, we obtain classical field 
equations: 

^1 + ^2 

v^v^x + (v^xtv'^x)x + ^^(v^V'tv'^x)^ = 0. (5) 

(Compare equations (63) of Q.) 



For cosmological purposes, we are interested in evolving very inhomogeneous 
field configurations, in regions where a texture is collapsing and the nonlinear 
terms in the equations are very important. There is little hope of making much 
progress analytically in general, and recourse to numerical techniques is needed. 

We wish to discretize the field equations on a grid. However, any direct 
discretization of the field equations (^ will fail immediately and drastically: 
adding small corrections Ai/), Ax according to a discrete version of (|^) in general 
violates the sigma model constraints \^p\'^ = \x\'^ = 1, tp'^X = 0. With those 
gone, the fields will no longer represent a member of SU(3), and the simulation 
becomes meaningless. 

Interestingly, the standard method of dealing with this problem: to discretize 
the action and enforce the necessary constraints with a discrete set of Lagrange 
multipliers, does not work here. This is a pity, because this technique works 
beautifully for SU(2) (see Pen et al, , eqns (67) and (68)) and yields a very fast, 
efficient algorithm for evolving the SU(2) field forward. This method fails here 
because one must simultaneously solve for four Lagrange multipliers, (A^,Aj^, 
and the real and imaginary components of A^^), at each spacetime point. One 
needs to solve four coupled nonlinear (quadratic) equations for four unknowns, 
at each point on the grid, at each timestep. (For SU(2), one only needs to solve 
single a quadratic equation.) This is numerically costly and complex - it must 
be done by iterative methods, and one has to worry about choosing the right 
one of the sixteen roots. 

In order to solve the SU(3) field equations numerically, it is helpful to re- 
state the problem in a form which more obviously respects the structure of the 
group. As written in (ph , the field equations describe the motion of two three 
dimensional complex scalar fields. Thus there are twelve real degrees of freedom, 
which are four too many for describing a member of SU(3). We have to try to 
control the remaining four dimensions by imposing constraints in some clever 
fashion. It would be far better to work with a formulation of the field equations 
which only included the eight real degrees of freedom that exist in the problem. 

Note that this problem of extra degrees of freedom has nothing to do with 
this being a field equation: the same problem would exist if we had only a 
single pair of complex 3-vectors, which were constrained to be of unit length 
and always perpendicular to one another. In fact, this is analogous to a very 
simple classical mechanical system. If the two vectors were real, instead of 
complex valued, then the single point system would have the same dynamics 
as a pair of iron bars, of unequal length, welded together perpendicular to one 
another, and with their point of contact held fixed, but otherwise free to rotate, 
as in figure n^. In this analogous system, the motion is trivial: there are three 
degrees of freedom, so the motion of this system is completely determined by 
conservation of angular momentum. 

The pair of bars conserves angular momentum because their Lagrangian is 
invariant under rotations in 3-space, SO (3). Going back to the SU(3) break- 
ing Higgs field, the Lagrangian is clearly unchanged by an analogous global 




Figure 1: An analogous simple mechanical system. 

SU(3) rotation, xp -^ Uxp, x ^^ Ux- This invariance instantly gives us eight 
conserved Noether charges, components of angular momentum on this internal 
SU(3) space. Since our fields only have eight real degrees of freedom, their dy- 
namics are entirely determined by the conservation of these Noether currents. 
Further, if we formulate the problem entirely in terms of the local member of 
SU(3), U{x), and the local internal angular momentum densities, then it will 
be written entirely in terms of the eight real field degrees of freedom, and the 
corresponding velocities, with no extra degrees of freedom to be dealt with. 
In terms of Noether currents, then, the field equations are 










(6) 



where the Noether currents JfJ are given by 

j; = Im [vl'^^T'^d^^l, + vlx^T'^d^x] ■ (7) 

Here the T° are the generators of SU(3); for sake of convenience we will always 
use a basis where Tr[T°T''] = 5°'^ . Here and below, we will use the convention 
that if X° are the components of a Lie algebra member, then X — X'^T"' G 
L(SU{3)) is the algebra element itself. 

Now, define the local internal angular momentum densities, L°- — Jq . In a 
flat FRW universe, the field evolution equations, (||), become 



drL" 



a 



Im {w2^tyay2^ ^ ^2^tr,^2^| _ 2«2,« . 



(8) 



Here g^^ — a^(r)diag{l, — 1, — 1, — 1}, so a{T) is the usual scale factor; we are 
working in conformal time and comoving spatial coordinates: we assume the 
background spacetime is flat Friedmann-Robertson- Walker. V^ is the ordinary 
flat space Laplacian. The lTii{di'4>* T°' diip} terms have vanished from (pi) because 
the hermiticity of T" makes the quantity in brackets purely real. 

With equation (||), we have a straightforward, easily calculable method for 
evolving the angular momenta, L, forward provided we have the fields xj} and 
X- In order to have a complete method of evolving the fields forward, we need 
to be able to evolve the SU(3) field forward, provided we know L. To go back 
to the mechanical analogy: we need to solve for angular velocity in terms of 
position and angular momentum; we need a moment of inertia tensor. Using 
the condensed notation of (||), we need drU . Now, since U{x) must always 
remain upon the group, we know that there exists some 0(x) G L(SU(3)) such 
that 

drU = i^U , (9) 

and we need to get 57 from L. The necessary transformation can be accomplished 
by taking the definition of L: 



L^ = J^ = Im {vlxP^T^ ^ + vlx^T^ x) , 
and constructing the quantities: 

-0tL'^r>, xP^L^T^x, xp''L''T''{xp'< X x'') , 



etc. 



(10) 



(11) 



These components may be evaluated, to solve for xp,x^ by using the following 
general relation: For any complex 3-vectors, A, B, C, 



Im [A^'T'^B] T^'C ^ U (B^'C) A - (A^C) B + \ (sU - A^B) C 



(12) 



a kind of complex generalization of the BAC-CAB rule for double cross products. 
Equation (|2) may be readily proven by writing A^T'^B = Tr [T'^B^'I'j ^ and 
then using the various known properties of the T°s. 

Using ([l2|), then, we can evaluate the nine quantities of ([ll|), and solve for 
fl. If we temporarily rotate coordinates so that [/ = 1, then if we write 
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L* 
L* 



'aip 



T * 
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(where a is a shorthand for ij}^ x x^)^ then Vt is 
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(13) 



(14) 



This can more be written more concisely as £7" = I\L'', where 1°^ is the (almost 
diagonal) moment of inertia tensor, whose action is defined by (H). 

Finally, (llj) is only valid when U = 1. In order to obtain the general relation, 
we must do and undo a rotation to the U — I frame: 

n{L, U) = U{lo (U^LU) } C/t . (15) 

And that completes the circle. Taking (|lq) , (|l4]) , (g), and (||) together, we 
have the field equations in a form where we have explicitly solved for the time 
derivatives of two fields, with no extra degrees of freedom. The equations upon 
which we will base the numerical field evolution are: 

if = in{L,u)u , (16) 

L'' = Im{t.2,/,tT'^vV + t^2X^T^''V2x}-2-L" , (17) 

where the precise form oin{L, U) is given by (15),([l4|). 



3 The numerical algorithm 

Now that we have the field equations, ( [iq ) and ([l7|), we can set about construct- 
ing a field evolution algorithm based upon them. Space and time are discretized 
into grids. Since U and L are a field and its time variation, it is natural to place 
them on alternate half timesteps. We need to construct a finite difference tech- 
nique for evolving each field forward, if the other is known one half timestep in 
the future. 

3.1 Discrete evolution of L 



One can read the algorithm for evolving L forward directly from (17): 
a'(7-„-i/2) „ 



ja 



+ 1/2 



O- [T^n+l/2) 



^^rimlviipiT'^ [^:t' + <-' + - - et^j^^"^*^) 



where the ... signify similar sums over nearest neighbors in the j and k directions. 
Here and elsewhere, if spatial indices are omitted, the (i, j, fc)"^ point is implied. 



The spatial derivatives piece of ( |18| ) is clearly second order accurate, and the 
factor I "^^""'"^''n J takes the effect of the expansion of the universe into account: 
angular momentum density generically redshifts away with a^^ as the universe 



expands. (That looks a little odd because we are in conformal time. In ordinary 
time, the angular momentum density would decay as a^"^, and integration over 
space with the proper volume element would give you a genuinely constant 
Noether charge.) Also, ( p^ has the pleasant property that it exactly conserves 
the total amount of each component of global angular momentum on the grid, 
since the sum over the entire grid of all of the discrete second spatial derivatives 
will be exactly zero (periodic boundary conditions are enforced). 

And so we have eight Noether charges which are conserved in the continuum, 
and which are exactly conserved by this scheme of discretely evolving L forward. 
Since this is only an exact global conservation, it says little about the accuracy 
of this method. (The local conservation of the internal angular momenta will 
only be correct to second order in A^, in general.) However, it does push 
heavily towards numerical stability, for it is more difficult for solutions of finite 
difference equations to diverge if they are required to separately conserve eight 
independent quantities while doing so. 

3.2 Evolution equations for U 

And so the evolution of the angular momentum field is well in hand. Next, 
we turn to the evolution of the SU(3) field, U . We would like to have a dis- 
crete version of the continuum equation, ([l6|). The first requirement is that the 
evolution preserve the fact that U e SU(3), and so we want something like: 

Un+i = e^^-""+i/^f/„ , (19) 

which would be second order accurate in time, and would preserve SU(3) mem- 
bership up to machine accuracy. However, in calculating (119), one immediately 
runs into a snag: the moment of inertia transformation from L ^ fl explicitly 
depends upon U. Thus, in order to calculate ^n+i/2 we need a value for U one 
half timestep into the future. So we need some kind of predictive method to get 
a value for C/„+i/2, at least first order accurate in A^. 

Now, one can immediately think of several ways to construct such a thing: 



for example, one could approximate CX,i+i/2 = \/f/nt^„_i Un- However, in this 
construction, as in all the rest, one is immediately confronted with ambiguities: 



is that last construction better, or worse, than taking U.n+1/2 — \/ f^„_it^n f^n? 
Or should one try some kind of mean of the two? And there are a myriad of 
other possibilities, using the exponentials of combinations of i„-i/2, in+i/2) 
together with field values. These all have similar ordering ambiguities. The 
central problem is that constructions of the form (^^ are trying to use group 
multiplication, instead of addition, to advance the fields. This is lovely, and 
keeps one on the group, but unfortunately in SU(3) group multiplication doesn't 
commute. There will inevitably be ambiguities about the order in which the 
multiplication is to be done. 



Instead of working on a group, and using multiplication, we should really 
work on a space in which addition is legal, and so no ambiguities of ordering can 
arise. The obvious such space is the Lie algebra: we should take the logarithm 
of the continuum equation dlq), and discretize that instead. If we write 



Uix) 



= piFix) 



Fix) e L(SU(3)), 



then (|lq ) becomes: 



d^e'^ =tn{L, e'^) e'^ 



(20) 



(21) 



Now we need to solve (^) for F- This may be accomplished by means of the 
following construction (well known in lattice gauge theory): 



-'^m e'^ 






f [-iFe-'^^dre'^^ + ie'^^dr{Fe'^^)] dA 
Jo 






dA 



-iAAd[F] 



o F dA by the definition of adjoint. 



iAAd[F] ^y^ 



OF 



'g-iAd[F] _ 2' 
I I -^^tt:^ I o F 



-ikd [F] 



F 



iM [F] 



g-jAd[F] _ I 

iM [F] 



-Ad[F] ^ ^ 



g.Ad[F] _ I 



on . 



(22) 



And so we have an explicit expression for F in terms of Q. Now, evaluation 
of (p2) is little tricky, since the 8x8 matrix e*'^'*'^] — 1 is singular, and (E2|) is 
written in terms of its inverse. There is no real problem, since the Ad [F] on 
top is also singular, and the two singularities cancel one another out. However, 
one cannot calculate the right hand side of (E2h directly; one must interpret 

( e'Ad[Fi_i ) as the appropriate power series in i Ad [F] . So long as the series 

converges, there is no ambiguity as to the meaning of (E2). 



We therefore have an expression for F, which can be expressed in terms 
of a power series we can find|j, of quantities that we know. At this point we 
could use a discrete version of (p^), expressed as a power series, and attempt 
to construct a field evolution algorithm based upon that. That would be a 
sound way to proceed, but slow, as the power series may take many terms to 
converge. (Indeed, sometimes it doesn't converge, as we discuss in appendix 
A.) Instead, with a little bit more maneuvering, we can construct an exactly 
calculable expression for the right hand side of (|2^). 

We proceed as follows: Write out the power series expansion of (|2S 

1 + i?i*Ad [F] + 1^ (zAd [F]f + 1^ {iM [F]f + 



o n (23) 



= f7 + Bl^[F,f7] + ^[F,[F,17]] + ^[F,[F,[F,f]]]] + ... , (24) 

by the definition of adjoint. Here the -B„ are the Bernoulli numbers. Now, 
diagonalize F, 

F = RFdR^^ , and define Q' = R^^flR . (25) 

With these terms, the power series (EJ) becomes: 

F = Rn'R-^ + Bii[RFDR^\Rn'R-^] 

■2 o 

+^ [RFDR-\[RFDR-\Rn'R-^]] + ... 



rIq' + iB,[FD,n'] + ^ [Fd, [FD,n']] + ...|i?-i 



(26) 



Now, we can express fl' in terms of complexified Lie algebra elements (step 
operators): 

n' = n'jj + {n'TT+ + %u+ + n'yV+ + h.c.) . (27) 

(The notation used here is standard for SU(3), see e.g. ref. [|| pp. 97-102. Also, 
we abbreviate $7'^ = Uip^T^ + flyY is the diagonal part of O'. The components 
of ^'jj are purely real; the rest are complex.) Now, in the complexified basis of 
the Lie algebra, the commutator of a diagonal matrix with any basis element 
is proportional to that basis clement. Thus we have a chain of identities of the 
form 

[ Fd,[Fd,...[Fd ,T+]]...] = a^^r+, ^28) 

n terms 

and similarly for all of the other basis elements. Here aj,^{FD) is a simple real 
valued function of the diagonal components of Fd. Using these relations, we 

•^Actually, the power series is very well known, since J'_ is a generating function for the 
Bernoulli numbers. 
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can write the power series as 



similar terms in other basis elements >R ^ . (29) 



And now we can explicitly resum the power series, which is now just a power 
series in an imaginary number. This gives 



«"[/ + 



r2'yt/+ + h.c.s \r-^ . (30) 



e'""+ - 1 

The diagonal components of 17' are uneffectcd because they commute with the 
diagonal matrix F^, so 0:^3 = a^ ={). All other components of the complexified 
basis of the Lie algebra have some nontrivial commutator with the Fd^ and so 
pick up an extra factor of i^ by this transformation. However, the a s are 
just real numbers, and so now this function is very easy and fast to calculate. 

And so we are done: in ( pO| ) (together with definitions in ([25|),(p|)) F is 
expressed entirely in terms of elementary functions of i^, \J = e*^, and L. 
The only fancy operation necessary is the diagonalization of a 3 x 3 hermitian 
matrix. Everything else is addition, multiplication of 3 x 3 complex matricies, 
and exponentiation and arithmetic of complex numbers. 

Everything in (^0|) can be calculated explicitly, and to machine accuracy, 
without resorting to iteration. And, since F has by construction precisely the 
eight physical degrees of freedom that exist in the problem, evolution of F is 
guaranteed to keep you exactly on the group. Finally, since _F is a member of 
an algebra, addition is legal, and so we may in fact justly discretize ( |30|) as a 
difference equation. 

There is only one more sensitive point which has not yet been addressed: 
the Lie algebra has a different topology than SU(3). Thus, one should be quite 
suspicious of any method of studying topological defects in SU(3) which looks 
at the group through the Lie algebra; the exponential map between the algebra 
and the group is not one to one. This point is explored in appendix A; the result 
is that one must be a little careful to always apply equation (BG) in a context 
in which it is meaningful. However, this does not materially change change 
the field evolution scheme. Also, note that the evolution equation for L, (ttq), 
depends only upon the SU(3) group members, U {— [..,i/?,x]), and not upon 
the Lie algebra member F . So there will be no problems with the evolution of 
the L field coming from the topology of the Lie algebra attached to F . 
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3.3 The finite diff'erence algorithm 

And so, we are finally in a position to construct a finite difference algorithm 
for evolving the SU(3) nonlinear sigma equations. The differencing scheme 
we use for L has been given; for F we need to use some scheme to predict 
F one half timestep into the future, since F depends upon F. To make this 
prediction accurately, we use a "— i,+|,+l" differencing scheme as follows. If 
we abbreviate (jsy) as 

F^G{F,U,L), (31) 

then, assuming Fn, Fn^i, Ln-1/2, Ln+1/2 a-^e known, we construct: 
3 1 

3 

Fn+l/2 = -Fn-1/4 + T^rG(i^„, [/„, 2(-^ri+l/2 + -^n-1/2)) (33) 

Fn+1 = Fn + ArG [Fn+l/2, Un+l/2,Lnj^i/2) ■ (34) 

And, again, the L field is evolved by 



_ g (^n-1/2) Ta 
1+1/2 a^(r„+i/2) "-1/2 






This difference scheme leads to an approximation to -F„+i/2 which is second 
order accurate in A,- if L is changing slowly. The difference scheme is shown 



geometrically in figure 3_^; two second order, the correct path between the 
points between the points {Fn^i, Fn, Fn+1/2) is a circle whose tangent is given 

by Fn- If angular velocity is changing reasonably slowly, the angle between 
Fn and i^n+1/2 should be half the angle between Fn and Fn-i, and a simple 
geometrical construction shows that a line parallel to the tangent at Fn will 
cross Fn+i/2 if started from the point labeled F„_i/4, as defined in ([3^). In the 
worst case, when the angular velocity about the circle is changing rapidly, this 
will be at worst 1^' order accurate in A,-, and so Fn+i will always be at least 
second order accurate in At-. 

Experimentally this works much better than any of the (faster) first order 
accurate schemes we used for constructing Fn+i/2- This scheme is somewhat 
costly in computation time: it involves two calculations of the complicated 
function G{F, U, L) per timestep, on each grid point. However, this differencing 
scheme, combined with the discrete L evolution equation (nsh (together the 
technicality discussed in appendix A), yields an evolution algorithm which is 
extremely stable (solutions do not diverge even when totally random initial 
conditions are given, in Minkowski space. Minkowski space is the hardest case, 
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Figure 2: The differencing scheme. To second order, we may approximate the 
correct path as a circle defined by F„, Fn-i, and the tangent Hne F„. The hne 
between -F„_i/4 and i^„+i/2 is paraUel to the tangent at Fn. 

since there is no damping coming from the expansion of the universe.) Further, 
it's behavior is very accurate by every test we have put to it. 

4 Testing the numerical code 

In the last section, we presented a finite differencing algorithm, equations (ttq), 
(p2f)-(p3), designed to approximate solutions to the field equations for a broken 
global SU(3) symmetry, which were developed in section pi We have subjected 
this algorithm to numerous tests of its accuracy and stability. This section 
summarizes these tests, and their outcomes. The central conclusion is that the 
algorithm is stable, under any set of initial conditions with ^initial — 0, provided 
the timestcp ^ < 0.3. Further, when the algorithm is stable, its accuracy (as 

measured by violation of local conservation of energy) goes as ( ^ ) . Also, 

when fed initial conditions where the analytic solution is known, or where the 
field behavior may be compared with the results of other programs, this program 
produces results consistent with the correct or previously found ones. Unless 
otherwise specified, all tests were done in Minkowski space, with a = 1, d = 0. 

4.1 Stability 

Firstly: as to stability. This algorithm has a strong tendency to be stable, 
since by the way it was constructed it will exactly conserve the eight global 
components of internal angular momentum. This does not, of course enforce 
stability, for this still allows a momentum component La — *■ cxa in one region and 
La —>■ — cxD in another region, the energy of the system diverging while preserving 
the correct total momenta. Experimentally, we find that the total energy of the 
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system tends to diverge roughly linearly in time when -^ > 0.4, and that total 
energy remains bounded, doing some random walk about a mean, for ^ < 0.3. 
Suprisingly, the simulation remains stable even when we are totally unfair to 
it, putting in initial conditions which are random group elements at each grid 
point, with no damping coming from the expansion of the universe. (In this case 
one cannot really interpret the system as an approximation to a field; rather it 
is an array of randomly oriented SU(3) members connected by springs.) All of 
these results come from runs where the initial momenta L are zero. If one puts 
in initial momentum, together with an initial field configuration consistent with 
that momentum, the system behaves predictably: its energy is stable unless one 
makes the initial momenta L so large that F A,- > 0.4 on a fair fraction of the 
lattice. (Putting in an initial field configuration Fo,Fi which is not consistent 
with the initial momenta yields very unstable behavior, as might be expected). 

4.2 Local energy conservation 

Nextly, as to accuracy. The best generic method we know of to test the accuracy 
of solutions is to watch the violations of local conservation of energy. (This is 
a much stronger test than merely looking at global energy conservation, where 
local errors tend to be washed away by averaging.) This check is a little bit 
tricky to implement, for to use it one must construct discrete derivatives of the 
discrete energy momentum tensor 9'^To^ which center properly upon a single 
grid point, and then check the extent to which 

^o^oo - SiTot = e (35) 



differs from zero. The particular discretization chosen for the components of ( p5D 
are written out in appendix B. We find that, with a timcstep of At-/A^ — 0.1, 
with fairly smooth initial conditions, the rms violation of energy conservation is 

M - -02 . (36) 

where the average is over grid points. So for smooth fields, average violation of 
local energy conservation is about 2% per grid point, for a timestep of .1 grid 
units. With totally random initial conditions, this went up to about 4-5% per 
grid point, per timestep. We found that this error scaled roughly with A^, as 
it should, and that it did not systematically tend in one direction or another. 
Thus on large grids, the global energy conservation looks ridiculously good, as 
many small errors average away. We found that for a fixed length of evolution 
time, the error in global energy conservation went with A^, for every class of 
initial conditions we tried. These results remain unchanged when we put in the 
expansion of the universe and added the necessary extra term to the the energy 
conservation equation. 
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4.3 Comparison with known behavior 

Next, we checked the behavior of the simulation in certain known circumstances. 
Firstly, we set the U equal to a constant over space, gave it a constant angular 
momentum density L. Not suprisingly, the entire U field rolls peacefully around 
SU(3) in unison; this was almost the only thing that it could do, since the 
angular momentum field L will never change under these circumstances. For 
a less trivial test, we set % = (1,0,0) everywhere on the grid, and gave the 
ip field some initial conditions ip = (0, a,/3). In this case a perfect evolution 
would never excite %; all motion would remain on the SU(2) subgroup defined 
by fixing x- This is also one circumstance in which results from this code can be 
compared with results from another code, namely that for SU(2), based on very 
different principles [^ . We found, firstly, that x is in fact very little disturbed: 
when we put in a single large texture on the ip field, and let it collapse, the x 
field was excited to about 10""^, on average, away from its starting point (by a 
time of 8 grid units). During the same period, the tp field ranged over its full 
extent, from —1 to +1 on the two components free to move, since the texture 
had collapsed by this time. 

Finally, we constructed the same texture in the tp field of this code (again 
leaving x fixed, so the field should be purely on an SU(2) subgroup) and on 
the previous SU(2) code. We evolved the two codes forward, with the same 
timestep, on the same sized grid; if both codes evolved perfectly, the single 
texture would collapse identically on both grids. The results of the two runs are 
quite dramatic: the two codes behave almost identically. 

The case of worst disagreement between the two codes happens just as the 
texture unwinds, when spatial gradients are largest and the two different dis- 
cretization schemes must show themselves against the grid. This worst case 
difference is shown in figure H, which are plots of the kinetic energy of the fields, 
on a two dimensional slice through the center of the simulation box. (The grids 
were 32 points on a side, and the texture was placed in the center of the box, 
with no initial kinetic energy. The fields were evolved forward with a timestep 
of .1 grid units.) The energy scales in the two pictures are the same; all num- 
bers are in grid units. Note that this discrepancy is only to be expected: both 
simulations will necessarily be inaccurate around the center of a collapsing tex- 
ture in the final instants of its collapse. The instant of texture collapse is, by 
definition, a time when neighboring grid points at the center of the texture lie 
upon extremely different places in the group. No finite difference approximation 
to a partial differential equation will accurately describe such a situation. 

Once the textures have collapsed and the gradients in the fields are smaller, 
both codes may be expected to be much more accurate. Indeed, later in the 
simulation, as the energy of the textures is radiating outward in a spherical shell, 
the two simulations are virtually identical, as shown in figure y. Thus when its 
initial conditions are restricted to an SU(2) subgroup of SU(3), the behavior of 
this SU(3) code is in excellent agreement with the SU(2) nonlinear sigma model 
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Figure 3: The kinetic energy density on a two dimensional slice through the 
center of a single texture, which was started with the same initial conditions on 
the SU(3) and SU(2) codes. This plot comes at the instant of texture collapse, 
when the field evolution will be least accurate. 



code which has been previously developed. 

5 Results 

This code was developed in order to study the effects of cosmic texture, coming 
from a broken global SU(3) symmetry, upon the early universe. In particular, 
we would like to use this code to predict the appearance of the cosmic mi- 
crowave background, if the primordial density fluctuations were seeded by this 
type of texture. In order to make predictions about the CMBR and density 
fluctuations, we need to evolve this fleld forward through the early universe, 
calculating parts of its energy-momentum tensor as we go: on the large scales 
of interest, this SU(3) fleld will interact only gravitationally with the matter 
and radiation of the early universe. After calculating the necessary parts of the 
energy-momentum tensor of this fleld, we can then use this information to pre- 
dict the power spectrum, nongaussianity, etc. of the CMBR, using previously 
developed techniques. 

In this section we describe the method of uniting this SU(3) evolution code 
with other codes which evolve the matter and radiation of the early universe, 
and give some results from these simulations. 

5.1 Scaling behavior 

In order to use this code to simulate the dynamics of an SU(3) field in the early 
universe, we need to construct a set of initial conditions. After the initial quench. 
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Figure 4: Kinetic energy density on the same slice of the SU(3) and SU(2) 
simulation, after the texture has collapsed. 

a physical SU(3) Higgs field would rapidly approach some scaling solution, first 
in the radiation era and then in the matter era. We do not know, a priori, what 
the power spectrum of the scaling solution would look like, and so we cannot put 
in realistic initial conditions. However, if we put in random initial conditions, 
and evolve the field forward, the simulation itself should also approach scaling 
behavior. So the first thing we must do is make sure that the SU(3) code does 
in fact approach a scaling behavior. 

In this case, we will look at the scaling behavior of one very important 
quantity for the growth of gravitational perturbations, namely the source for 
Newtonian gravity, p + 3P ex Too + Tu. This acts as the gravitational source 
term in the scalar density perturbation equations in synchronous gauge (See e.g. 
01 Hi)- This term alone is sufficient to determine the 'intrinsic' perturbations 
to the CMBR (the synchronous gauge surface terms) and more generally, the 
matter perturbation power spectrum. A brief calculation shows that for this 
SU(3) system, 

,2^ 



Too +Tu^2{ vfi> 



. 2 



V2X 



By scaling, we mean the scaling of the power spectrum of this quantity, 

1 



{\Mr)n 



where 61. — 
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{x,t)e' 



-ik-x j3 



d-'x 



(37) 



(38) 



is the fourier transform of (p. If the simulation is purely in the radiation or 
matter era, there are no length scales in the system (save the unphysical box 
size and grid size), so this power spectrum should reflect the scale free nature 
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of the problem. Thus, in the radiation era, we require 

(l0.(r)P) = ^ (39) 

for some function /{kr). (The extra factor of t is necessary to balance the 

. 2 
units; 4>{x) ^ tp has units of r^^, and there is a factor of t^ coming from the 

expansion of the universe). 

Numerically, we take the average in (139) by averaging over the power spectra 
(pl gotten from an ensemble of runs. To test the scaling, behavior we did an 
ensemble of six runs on 96^ grids in the radiation era, starting from random 
initial conditions evenly distributed over SU(3). As can be seen from figure g, 
the power spectra in a very wide range of wavenumbers k rapidly approach a 
scaling solution, and remain on it after r ~ 6 grid units. 

Since scaling in this power spectrum of (f) begins at r ^^ 6, we can find the 
scaling function /(fcr) by averaging the values of t (|(/)fc(T)p), for all fc, for times 
greater than six grid units (and less than half the box size, r < 48 in this case. 
For times later than this, the box size introduces an unphysical time scale into 
the problem.) The result of this averaging is shown in |g, which was constructed 
by summing all the values of r (|(/)/j(t)P) into bins in kr. 

And with this scaling function, /(kr), constructed, we are finally in a po- 
sition to begin comparing the behavior of SU(3) textures in the early universe 
with their SU(2) counterparts. A similar set of six runs, on 96'^ grids, in the 
radiation era, was done using the SU(2) code. Using the same analysis, we con- 
structed a scaling function fsu{2)i^''')- ^he two functions are shown together 
in figure ul The resemblance is striking: up to a multiplicative constant, the 
two functions are identical to within one another's error bars. The multiplica- 
tive constant is immaterial; it would be removed by normalizing both CMBR 
predictions to COBE, for example. Thus, at least as far as one can tell from 
looking at power spectra in the radiation era, SU(3) based textures have almost 
exactly the same effect upon the early universe as SU(2) based textures do. 

5.2 CMBR predictions 

Finally, we used this code to provide the calculate the appearance of the cosmic 
microwave background, in a universe where the large scale structure was seeded 
by SU(3) textures. In the simulations, we start out the field with random initial 
conditions, and let the field evolve forward in the radiation era long enough for 
scaling behaviour to set in. We then let the universe transfer from radiation to 
matter dominated epochs, and continue forward in the matter epoch until the 
time of last scattering. During this simulation, we used the energy-momentum 
tensor from the SU(3) field as a gravitational source for the density and velocity 
fluctuations, using code developed by Crittenden and one of us {m, M). 

From the density and velocity perturbations in the primordial gas of one run, 
in a 256'^ box, we made a set of 48 independant maps of the visible temperature 
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2,10,20,30,40 




Figure 5: The power spectrum of <f) approaches scahng in the radiation era. 
This is a plot of r (|</'fc(T)p) vs. A:t, each scaled to remove the box size. Each 
trace in shows the time evolution of a different fourier mode; each starts at zero 
because the simulation is started with the fields static, so (/)(a;, 0) ~ ■0^ = 0. 
Each point on each trace was taken at a time r one grid unit later than the 
preceding point; thus each trace starts at r = 0, has its first joint at r = 1, and 
so on. (Since the grid is 96 points on a side, the fourier wavenumbers run from 
k = (1-48)-^. Thus the range of k shown covers almost all the wavelengths 
present in the simulation). Each of these fourier modes represents an average 
of the power spectra from six separate runs. In each run, the power spectrum 
for a given |fc| is an average of power spectra in three perpendicular directions. 
Thus each point here is an average of eighteen independent quantities. 
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SU(3), t = 6-40 




Figure 6: The averaged scaling function fikr) for SU(3) in the radiation epoch. 
This is a plot of /{kr) vs kr, and was constructed by averaging all values of 
T (|(/)a;(t)|^) into bins in kr. The average runs over all values of k, for t in a 
region where the power spectrum is scaling, 6 < r < 40. The vertical error bars 
come from assuming that (j)k{T),4'k'{T') will be independent from one another 
ii k ^ k' . Fourier components with the same wavenumber, at different times, 
are not assumed to be independent. The horizontal error bars are just the bin 
sizes. 



1 

(h-5- 


SD(2) (above) 


SD(3) (below 


, log-log 


-0.5 


0.5 




1.5 2 


-1.5 






X 


-2 






^% 



Figure 7: A comparison of the scaling functions fsu(2)(^''') ^^'^ fsu(3)i^''')- 
This is a plot of logj^Q / vs. logj^Q kr for both scaling functions; the SU(2) 
scaling function is the upper one. Both of these scaling functions were created 
by averaging power spectra from six radiation era runs on 96"^ boxes. 
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fluctuations on the surface of last scatter. Two such maps, which would appear 
as 10° X 10° patches of sky from earth, are shown in figure |l Superficially, the 
CMBR anisotropy maps presented here are very similar to previously published 
textures maps. (See, for example, the pictures in M.) Note that these maps only 
represent the intrinsic anisotropics generated on the the surface of last scatter; 
they do not include the effects from propagating the light from the surface of 
last scatter to us. This omission has little effect upon the small scale structure 
of the image, but it does mean that the images have too little power on large 
angular scales (l < 200). 

From this assembly of images, one can readily calculate the high end of the 
CMBR power spectrum, which is shown in figure ^ (along with similar plots 
from SU(2) textures, cosmic strings, and monopole calculations). Note that 
these maps don't contain the necessary large scale information to normalize 
to COBE, so the normalization of ^ is left arbitrary. We do not expect that 
the large scale behaviour of the CMBR spectra will be much different between 
the SU(2) and SU(3) models, when it is calculated: the interesting effects of 
the nonlinear dynamics show up on small scales. In particular, the positions 
and shapes of the doppler peaks are correct in this plot. Note that power 
spectum from the SU(3) textures has a high, second peak, where the SU(2) 
power spectrum has none, although the position of the first peak is the same. 
So the SU(3) field dynamics give a bit more power on smaller scales than does 
SU(2). 

Finally, note that the CMBR maps of contain a handful of extremely high 
peaks, many more than would be expected in a picture of gaussian fluctuations. 
These high peaks mark the infall of gas into the gravitational wells created by 
the energy of the collapse of individual textures: they are a primary source 
of primordial density fluctuation in this theory. The primary distinguishing 
features of the textures CMBR maps is an overabundance of large (~5 sigma) 
peaks, of which there are 0(1) in each map we created. To quantify this, we 
counted the number of peaks, of a given height, which appear in a given area 
of sky in the collection of maps that we created. The distribution of peaks for 
SU(3) and SU(2) textures is shown in ^. Somewhat suprisingly, this measure of 
the nongaussianity of the maps gives almost indistinguishable results for SU(2) 
and SU(3) textures. Thus the density of collapsing textures in the universe is 
almost identical for the two theories. 



6 Conclusion 

We have developed and used a reliable numerical technique for studying the 
effects of textures from a broken global SU(3) symmetry upon the early universe. 
We find that, on the whole, SU(3)-generated textures behave very similarly to 
textures which come from breaking a global SU(2) symmetry. In particular, the 
scaling behaviour of the two theories in the radiation era, and the density of 
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Figure 8: Temperature anisotropies on the surface of last scattering. The angu- 
lar scale of these maps is 10° on a side, for fl — I, flB = 0.05, and h = .5. The 
color scale units are in one standard deviation a for the CMBR anisotropy. 
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Figure 9: The region of the CMBR anisotropy power spectrum dominating 
on small angular scales for SU(3) textures (filled boxes), SU(2) textures (open 
boxes), and global monopoles (open triangles). All curves are for fi — 1, h — .5. 
Error bars show the statistical uncertainty (one sigma) in the SU(3) result. 
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Figure 10: Differential number of local maxima (and minima) in the CMB 
anisotropy, of height va (and —va) per square degree, where a is the rms. The 
number density of maxima is shown in solid lines, minima in dotted lines. The 
filled boxes show the SU(3) texture results, the open boxes the SU(2) results 
(from ref. [3]). 
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nongaussian peaks on the cosmic microwave background, are identical in the two 
theories to within our numerical errors. We found that the CMBR power spectra 
of SU(2) and SU(3) textures, are close, but do have a measurable difference 
in shape: the SU(3) CMBR power spectrum has a bump not present with 
SU(2) textures. The two textures theories are much closer to one another than 
they are to qualitatively different theories of large scale structure formation. 
This suggests that the cosmological predictions of texture theories are rather 
robust: changing the particular group which is broken to produce the textures 
has relatively minor cosmological effects. However with the great accuracy one 
now anticipates in measurements of the CMBR anisotropy power spectrum, even 
these small differences may be subject to experimental test. 
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A The problems from nonuniqueness of U ^ F 

In section pi we noted that there were some subtleties in the evolution of the 
field F, the logarithm of the SU(3) field U, springing from the non-one-to-one 
nature of the map between the group and its Lie algebra. In this appendix, 
we discuss precisely how these subtleties arise in the context of this numerical 
evolution, and how they are handled. We show that, if treated properly, this 
ambiguity never presents a problem in the field evolution or interpretation of 
results. 

Recall that the continuum equation for the time derivative of F was of the 
form: 

F^ R Ui'd + /c!^l^_ ^ ^'t+T+ + etc. I i?-i , (|0|) 

where Vt' was some Lie algebra element, i7^ is its diagonal part, and the rest is 
split up in terms of the complexified basis. Here F has been diagonalizcd into 
RFdR~^, and the a^^, etc. are simple linear combinations of the components 
of Fjj. In particular, if we break up F^ into the usual generators of the Cartan 
subalgebra on the main diagonal, Fd — FtsT^ + FyY, then the as are: 



1 I 



Ft3 , a,j^ ^Fy - -Ft3 , av+ = Fy + -Fts ■ (40) 



For all the generators, a^_ = —a^^. 

Now, the expression for F, (pO|), is written out in terms of functions of the 
form -n; — 7 of the six a s. The behavior of this function needs to be watched 
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carefully, since our time derivatives are expressed in terms of it. It has a remov- 
able singularity at the a — 0, which causes no trouble: the function smoothly 
approaches i in that region. However, the function has genuine singularities at 
±2n7r, for all integers n other than zero! Thus if we naively evolve the F field 
forward according to (pW) we rapidly get into trouble: the time derivatives go 
infinite the first time one of the a s hits ±2tt. 

And so direct application of (|30|) only makes sense when we are within a 
certain neighborhood of the origin in the Ft3, Fy plane. The situation is il- 
lustrated in figure [ij: the region of good behavior of all a s is bounded by a 
hexagonal wall of singularities of -t^t — - in one or more a. 
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Figure 11: The singularities of equation (30), in the plane of the (T3, Y) compo- 
nents of Fu- The dark lines denote places where some a = 27r, so a(e'" — 1)~^ 
diverges. The dotted circles denote regions inside which a(e'" — 1)~^ < some 
chosen bound, for all a, for a particular choice of origin. 

Now these singularities are not physical, of course; the field equations written 
in the form of (Hq) , purely in terms of C/, are not pathological anywhere on the 
group. The singularities in (5fl) are just symptoms of the fact that, when F 
is far enough away from the origin, the coordinate system of the Lie algebra 
element F is not a very good one for this type of calculation. 

The solution to the problem is simply to change origins, whenever we get far 
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away from F ~ 0. That is, let us generalize our notation slightly, so that the 
physical field U is 

U ^ e'^O , (41) 

where O is some origin, which we may choose to alter from time to time as we 
advance the fields forward. The generic effect of adding an origin as in (^) is 
to change F and its diagonalization: 

Fo^Fd + ^d . R^R' . (42) 

We are permitted to change our origin in any way we please, and so we can al- 
ways choose the diagonal matrix A^i in such a way as to make Fd smaller. Pro- 
vided we monitor the values of Fd , and change origins whenever it approaches 
a singular region, we can always avoid any difficulty. 

In fact, figure [ll| suggests a particularly nice set of origins to use for this pur- 
pose. We would like to choose the new origins to be places which are as singular 
as possible in terms of the old coordinates, so as to keep track of the smallest 
possible set of origins. As such, the points A and B on figure ^ where two 
lines of singularities intersect, naturally suggest themselves. (The three points 
labeled A on figure [III actually correspond to the same point in SU(3), as do the 
three points labeled B. Every member of SU(3) can be made to correspond to 
a point in the interior or on the boundary of the hexagon. Within the hexagon, 
there is a three-fold "refiection" symmetry coming from the ambiguity of how 
we choose to order the eigenvalues of F within Fo)- In fact, the SU(3) mem- 
bers gotten by exponentiating A and B are just e^'"/^ , e^*^/^ times the identity: 
together with 1, they form the center of SU(3). These points are particularly 
nice origins because they make changing origin very easy; they commute with 
everything. In particular, the transformation of the R matrix in ( [42| ) becomes 
trivial, whereas in general it is a painful thing to compute. 

So: the evolution algorithm described in section || is modified as follows. 
Each time we diagonalize F, we construct the radius^ F^g -|- Fy and figure 
out how far we are from our current origin in figure O. If we are outside of 
a certain radius, we declare it time to change origins, and shift the origin of 
F to the corner of the hexagon nearest to the current position of Fd- (We 
must also shift the origin of the previous timestep in the same way, because the 
algorithm uses the previous timestep directly with this one.) The only thing 
that need be decided systematically is which cutoff circle to use. From figure 
[ll| we see that the radius must be greater than ^, so that the three circles 
centered at 1, A, B intersect, and less than 27r, to avoid the nearest singularity. 
We chose ^, which provides a nice region of intersection, but is far enough 
away from the nearest singularities that the time derivatives never become too 
big; (37r/2)|(e^*'^/^ — 1)~^| ~ 3.3 is the largest value obtained by the diverging 
function in that region. (The smallest value obtained is 1, at the origin). Now, 
since the entire calculation of F will always be within some circular region in 
which F is well behaved, we never encounter any problems. 
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Everything else proceeds as before. Of course, we must now keep track of 
which of the three origins associated with each grid point, and make use of the 
origin as in mV) whenever we construct U. Fortunately, the simple form of 
the three origins chosen means that this does not slow down the program very 
much. 



B Discretized derivatives of T^^,. 

In section ]q, we describe the constructing the discrete derivatives of the energy- 
momentum tensor, in order to test how well the numerical algorithm obeys local 
conservation of energy, d^T^o = 0. When discretized, this becomes 

SqTqo — SiToi — e , 



where e is the error we seek to measure. Now, (35) involves a time derivative of 
the 00 component of T. A discrete version of Tqo lives most naturally upon a grid 
point on a whole timestep. The momentum densities Toi involve constructs like 
dxip dotp, and so most naturally live on the links between grid points, though 

on whole timesteps. (We place the time derivatives, i/?,^ on whole timesteps 
since we need to use both U and L to construct these. The L grid lives on 
half timesteps, while U lives on whole timesteps, so one or the other will have 
to be averaged to produce a [/ at a single point. Since L lives in an algebra, 
averages of two L are meaningful, whereas naive averages of two values of U are 
not. Thus we chose to construct time derivatives which live on whole timesteps.) 
Since equation (Bq) itself involves a time derivative of Tqq and spatial derivatives 
of the Toi, the discrete version of it may be fairly easily centered on a whole 
grid point, on a half timestep. 

(43) 

(44) 

with a parallel term for x- Here the if) terms are calculated from the L field, 
using the mean between the two half steps, as in (p3). Next we need to con- 
struct a discrete version of d^T^x which lives on a whole grid point, on a half 
timestep. Since the Tq^ live on whole timesteps, on links between lattice sites, 
this must be constructing by differencing adjacent links, and averaging between 
two timesteps. In one direction, this gives: 



Tqq - 


= ^(T^or-^o'^) , where 


rpn 
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-|- similar terms in j, k directions) 
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Jox -Tq^ ]} ^ where (45) 



T^^'^' = ^-?Re U:, + ^;+M (C+^-V^:.) , (46) 



2A, 

with, again, a similar term for x- 

These are the explicit forms used for the tests of local energy conservation 
used in section |[ They are shown here in detail because there several possible 
ways of constructing these discrete T^i/ components and their discrete deriva- 
tives, and most of them do not work. 
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